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Abstract 



We calculate in (two-loop) chiral perturbation theory the local NN-potentials gen- 
erated by the three-pion exchange diagrams proportional to g\ and g\. Surprisingly, 
q ' we find that the total isoscalar central 37r-exchange potential vanishes identically. The 

individually largest 37r-exchange potentials are of isoscalar spin-spin, isovector central 
and isoscalar tensor type. For these potentials simple analytical expressions can be 
CN ' given. The strength of these dominant 37r-exchange potentials at r = 1.0 fm is 4.6 MeV, 

2.9 MeV and 1.4 MeV, respectively. Furthermore, we observe that the spin-spin and 
tensor potentials due to the diagrams proportional to g\ do not exist in the infinite 

J> ' nucleon mass limit. 

-XT ' 

in 
o 

CM 1 PACS: 12.20.Ds, 12.38.Bx, 12.39.Fe, 13.75.Cs. 

In a recent work |J we have started to calculate the static nucleon-nucleon potentials 
generated by the (two-loop) 37r-exchange diagrams with interaction vertices taken from the 
leading-order effective chiral 7riV-Lagrangian. In ref. ||] we restricted ourselves to the eval- 
^ \ uation of four relatively simple classes of diagrams which in fact comprise all graphs with 
the common prefactor g\/f®- It was stressed in ref. ]I] that the chiral 37riViV-contact vertex 
depends on the choice of the interpolating pion-field and therefore one has to consider rep- 
resentation invariant classes of diagrams (classes I and II in ref.[|l[]) by supplementing graphs 
involving the chiral 47r-vertex. Another obvious consequence of this is that there is no unique 
way (in chiral perturbation theory) of denning the so-called "pion- nucleon form factor" which 
is often introduced in phenomenological models. The only meaningful separation of the in- 
teraction potential is the one into point-like l7r-exchange and 3"7r-exchange. The four classes 
of chiral 37r-exchange diagrams evaluated in ref.|l| gave rise only to isovector spin-spin and 
tensor potentials. The strength of the resulting coordinate-space potentials was very weak 
with values of ±0.1 MeV and less at a internucleon distance of r = 1.0 fm. Compared to the 
chiral 27r-exchange potentials ||, || these are negligibly small corrections. 

The purpose of this work is complete the calculation of the chiral 37r-exchange NN- 
potentials and to present results for the diagrams carrying independent prefactors g\j f® 
and g\j '/®. We will see that a grouping of these diagrams into five different classes is very 
useful. 

Let us begin with some basic definitions in order to fix our notation. In the static limit 
and considering only irreducible diagrams the on-shell NN T-matrix takes the following form: 

T NN = V c (q) + V s (q) &i-ff 2 + V T (q) B x ■ q a 2 ■ q 



+ 



W c (q) + W s (q) ai-a 2 + W T (q) a x -qa 2 -q 
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where q — \q\ denotes the momentum transfer between the initial and final state nucleon. 
The subscripts C, S and T refer to the central, spin-spin and tensor components, each of 
which occurs in an isoscalar (V) and an isovector (W) version. As indicated, the (real) NN- 
amplitudes Vc(q), ■ ■ ■ , Wt{q) depend only on the momentum transfer q in the static limit. 
We are here interested only in the coordinate-space potentials generated by certain diagrams 
in which three pions are simultaneously exchanged between both nucleons. For this purpose 
it is sufficient to calculate the imaginary parts of the NN- amplitudes Vc,s,t(<i) and Wc,s,t(q) 
analytically continued to time-like momentum transfer q = ifi — + with fi > 3m,,-. These 
imaginary parts are then the mass-spectra entering a representation of the local coordinate 
space potentials in form of a continuous superposition of Yukawa functions, 

V c (r) = [°° dupe-rimVcfa) , (2) 

V s (r) = — / dfifie-^lrnVriiti) - 31mV s (i{J.) , (3) 

D7T r J3m w 1 J 

Mr) = tt^-, / rf/x/xe-^(3 + 3/xr + A 2 )Imy T (^) . (4) 

07T T J3m.Tr 

The isoscalar central, spin-spin and tensor potentials, denoted here by Vc,s,t{t), are as usual 
those ones which are accompanied by the operators 1, B\ ■ a 2 and 3<7i - r c?2 • f — B\ ■ ai, 
respectively. For the isovector potentials Wc,s,t{t) a completely analogous representation 
holds. 

Let us now turn to the evaluation of the NN-potentials from the two-loop 37r-exchange 
diagrams. Application of the Cutkosky cutting rules gives the relevant imaginary parts 
Im Vc,s,T{in) and lmWc,s,T(i^) as integrals of the squared NN — > 3n transition ampli- 
tudes over the Lorentz-invariant three-pion phase space. Some details about these techniques 
can be found in ref. 0] where the same method has been applied to calculate (two-loop) 
spectral-functions of nucleon form factors. The pertinent three-body space integrals are most 
conveniently performed in the three-pion center-of-mass frame. The corresponding on-mass- 
shell four-momenta of the three pions read in this frame: k\ = (ux,ki), k% = (^2,^2) and 
fcg = (/i — ui — UJ2, —ki — A^)- The mass-shell condition /c| = m\ determines the cosine of the 
angle between k\ and k 2 (called z) as 

zk x k 2 = uj x uj 2 - n(ux + u 2 ) + ^(/i 2 + ml) , fc lj2 = sju\ 2 - m\ . (5) 

In the chosen reference frame the four-momenta of the external nucleons ca n be taken as 
P\ = (A i /2?P , w) and = (—(J,/2,pv), where v is a (real) unit vector andp = yj '// 2 /4 — M 2 = 
iM + 0(M~ l ) in the heavy nucleon limit (M — > 00). Of course the assignment of and P 2 U 
to the external nucleon lines has to be done such that four-momentum conservation holds for 
the transition NN — >• 3n — >• NN. The possible nucleon propagators occuring in diagrams are 
always of the form 

2M 1 2M -1 

(6) 



( Pl _^.) 2 -M 2 it.kj-e {P 2 + k J f-M* tv-ki + e 



The sign of the e-terms in eq.(6) is uniquely fixed from the original relativistic expressions 
and in the heavy nucleon mass limit (M — > 00) the positive quantity e becomes infinitesimally 
small. With these techniques at hand it is possible to perform analytically at least two of the 
four non-trivial 37r-phase-space integrations ||]. As a check on the heavy nucleon formalism 
employed we recalculated the imaginary parts of the isoscalar electromagnetic and isovector 
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axial nucleoli form factors and found perfect agreement with the results of ref.Q. As a 
further check we applied the methods to calculate the imaginary parts of the 27r-exchange 
NN-amplitudes. In that case only a much simpler two-body phase space integral needs to 
be evaluated and one reproduces indeed exactly the results of ref.|2|, | in a rather short 
calculation. 




Fig.l: 3n-exchange diagrams of class V proportional to g A . Solid and dashed lines represent 
nucleons and pions, respectively. 

Let us now turn to the results for the (two-loop) chiral 37r-exchange diagrams proportional 
to g\ and g\. We start with the diagrams of class V shown in Fig. 1. As stressed in ref.0 
diagrams involving the chiral 37riViV-vertex or the chiral 47r-vertex depend on an arbitrary 
parameter a and therefore one should consider the full class V as one entity. Obviously, the 
last two pion-pole diagrams in Fig. 1 contribute via coupling constant renormalization also to 
the point-like l7r-exchange. This effect is however automatically taken care by working with 
the physical 7riViV-coupling constant g n N- From an inspection of the spin- and isospin factors 
occurring in the diagrams of class V one finds immediately that only non- vanishing isovector 
spin-spin and tensor NN-amplitudes Ws,t will be obtained. We find the following imaginary 
parts from class V, 

ImW^M = fjdu ld u 2 {kl + w + 3K - + |) } , (7) 

Z 2 <1 



Im W? ] 



(6/i 2 + 2mJ)(u 1 + w a 



m 



- ,u(V + 3m£ 



2<1 
2 i Q™2 N 



arccos 



(8) 



The inequality z 2 < 1 defines the kinematically allowed (Dalitz) region in the co>iu;2-plane 
(which is bounded by a cubic curve) together with the obvious kinematical constraints m„ < 
^1,2 < A* — 2> m n an d 2m,r < W\ + 0J2 < jtt — m n . Note that the same integrand as in eq.(7) for 
was found in ref. || for the spectral-function of the nucleon isovector axial form 
factor. In the chiral limit m n = one can evaluate the remaining double-integrals in eqs.(7,8) 
using the substitution w\ = fi(l — xy)/2, u 2 = ^y/2 which maps the unit-square < x, y < 1 
onto the (in the chiral limit) triangle-shaped Dalitz region. 

One finds from class V repulsive isovector spin-spin and tensor potentials with an r~ 7 - 
dependence, 



9A r 



16 



85\ 



(47r)5/. 6 V21 71 " + 36/ 



r) 



4 

9aT 



m^=0 



Wft V 72 



245 



(9) 
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Next, we come to the diagrams of class VI shown in Fig. 2. The isospin factor of the first 
and second graph is At\ • r 2 — 6 while that of the third and fourth graphs is — 4fi -72 — 6. 
The last two graphs are irreducible whereas the first two contain the iteration of l7r-exchange 
and 27r-exchange (triangle subgraphs). This iterative part has to be separated off in the 
construction of the NN-potential. A closer inspection reveals that the two types of diagram 
differ only by some signs in nucleon propagators. The first two diagrams carry a factor 
— [(iv-ki + e)(iv-k 1 — e)]~ 1 in comparison to a factor {iv-ki + e)~ 2 from the last two (irreducible) 
diagrams. The irreducible part of the first two diagrams is obtained by switching the sign of 
one e-term (— e — > +e) such that the expression agrees with that of the irreducible diagrams. 



Fig. 2: 3tt -exchange diagrams of class VI proportional to g\. Diagrams for which the role 
of both nucleons is interchanged are not shown. They lead to the same contribution to the 
NN-potential. The (irreducible) isoscalar NN-amplitudes sum up to zero. 



In order to make this procedure more understandable consider the following integrals: 
J\dx[(x + ie)(x-ie)]- 1 = vr/e - 2 + £>(e 2 ) and f\ dx(x + ie)~ 2 = -2 + 0(e 2 ). The difference 
between both diverges as 1/e for e — > + . According to the definition in eq.(6) one has 
e ~ M~ l and from the (non-relativistic) Lippmann-Schwinger equation H it is known that 
the iteration of the potential leads to a contribution proportional to the nucleon mass M. 
Therefore the 1/e-term which gets subtracted by switching the sign (of one e) corresponds 
indeed to the iterative part. In the case of 27r-exchange all this has been worked out in detail 
in ref . || and as already mentioned the present methods allow to reproduce exactly the results 
of ref.p] for the irreducible 27r-exchange. 

After subtracting the iterative part the first two and the last two diagrams in Fig. 2 become 
equal up to a minus-sign. Combining this with the isospin factors one obtains again only a 
non- vanishing contribution to the isovector spin-spin and tensor NN-amplitudes Ws,t- We 
find the following imaginary parts from class VI, 



lmW^ VI \i/2) 



■2g 



2\3 



du)\du) 2 \ — k\ 



5 , n \ / 

-fjuvi + (imjJi - raj [z + 



arccos 



2 2 <1 



(10) 



fi fJ, {07TJ 7T J J J I OK! 



f 2wi 



+ 



(/i 2 + m 2 )(/i - 2ujj 



2o; 2 ) + 4/icuiu;2 



Z 2 <1 

arccosf— z) 



k 2 

2zuj 2 h 3c<Ji — 2 ii 

ki 
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In the chiral limit (m n = 0) one gets now a repulsive isovector spin-spin potential and an 
attractive isovector tensor potential of the form, 



4 

9A r 



(4tt) 



-(—) 

fl V 36 J ' 



W. 



(VI), 



4 

9A r 



=0 (47r)5/. 6 V3 
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Next, we come the diagrams of class VII shown Fig. 3. The isospin factor of all four 
diagrams is 6 and after considering their spin-structure one finds immediately that there will 
be only a non-vanishing contribution to the isoscalar spin-spin and tensor NN-amplitudes 
Vs,t- In the case of class VII one can actually solve all integrals analytically and one obtains 
the following closed form expressions for the imaginary parts, 



- ^ 



— V 



Fig. 3: 3ir -exchange diagrams of class VII proportional to g\. The isospin factor of these 
diagrams is 6. 



ImV, 



lmV^ VII \tfM) 



(VII) 



g\(fx - 3m,, 



357r(32/3)J 



2ml ~ 12^m w - 2^ 2 + 15^ + 2^f + 3^f 



g\(fx - 3m*) 



- 2,nn: I 6m* • In— 9^f 27^f 



(13) 
(14) 



,.3 , q,.2_ , o..™2 , 

35^(32^3)2 ' 

It is even more astonishing that the corresponding coordinate space potentials (inserting 
eqs. (13,14) into eqs.(3,4)) can be expressed through a simple exponential-function multiplied 
by a polynomial. We find the following repulsive isoscalar spin-spin and tensor potentials 
from class VII, 



yf"> W 



4 — 3?T?,7r V 

'' 1 — (1 + m*r) 2 (2 + m n r) 2 , (15) 



2(8tt/2)s r' 
g\ e - 3m ^ r 
2(8tt/2)3 ' 



1 + m n rY(l + m w r + m 2 r 2 ) . 



(16) 



f - -x- - > 



Fig. 4: Sir- exchange diagrams of class VIII proportional to g\. 

Next, we consider the diagrams of class VIII shown in Fig. 4. The isospin factor of the 
first graph is 7fi • f 2 — 6 while that of the second one is 1t\ ■ r 2 + 6. In order to separate 
off the iterative parts from the first diagram one has to switch the sign of two different e- 
terms in nucleon propagators. After this procedure the factors coming from the nucleon 
propagators agree identically for both diagrams in Fig. 4. The two diagrams differ however 
in the ordering of the first and third (spin-dependent) pion-coupling to one nucleon line. 
Exploiting furthermore the properties of a-matrices one finds that the spin-independent part 
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of this particular 37r-exchange interaction is equal with opposite sign for both diagrams, 
whereas the spin-dependent part is equal with the same sign. Combining this result with the 
isospin factors one obtains only non-vanishing isoscalar central (Vc) and isovector spin-spin 
and tensor (Ws,t) NN-amplitudes. Explicit evaluation of the diagrams of class VIII leads to 
the following imaginary part of the isoscalar central NN-amplitude, 



'c 



M = 'W^kf [h' du A l + 7T=f axccos(-z)} . (17) 



In the chiral limit (m n = 0) one finds a repulsive isoscalar central potential with an r 7 - 
dependence of the following form, 

% VnI) {r) = 7^4(60 - 4tt 2 ) . (18) 

A similar calculation gives for the imaginary part of the isovector spin-spin NN-amplitude W$ 
generated by class VIII the following result, 



ImWf^/i) = J 0^JJdu l du 2 [\{ul-^-Aml + ^u Jl ) 

Z 2 <1 



+ 



Z + 



arccos(— z) 



2 



[lux - m% 



2 r -v ™2 



2 m; - (jw 2 



(19) 



fci(l -z 2 ) 

However, as indicated by the question mark this results is of no use. On the boundary of the 
Dalitz region, z 2 = 1, the integrand in eq.(19) has singularities of the form (1 — z 2 )' 1 and 
(1 — z 2 )~^l 2 such that the double-integral JJ z2<1 dcuidcu2 diverges. The same non-integrable 

singularities were found in an explicit calculation of ImW^ m '(i^). Therefore one has to 
conclude that the isovector spin-spin and tensor potentials generated by the diagrams of 
class VIII do not exist in the infinite nucleon mass (M — > oo) or static limit. Note that this 
statement applies separately to the second (irreducible) diagram in Fig. 4 and the occurrence 
of these singularities has nothing to do with adding the (irreducible components of the) first 
diagram. Presently, we do not have a detailed understanding of the problem encountered here. 
Most likely it originates from the fact that several nucleon propagators in the infinite mass 
limit M — > oo can introduce severe singularities in loop-functions as the following one-loop 
example shows, 

dH l l 



/ 



i(2vr) 4 (l - uf{ml - I 2 ) 8tt 2 



m —lu 
+ ~, — n aN0/o arccos 



(20) 



m^ — uj 2 {m 2 — u 2 ) 3 / 2 m n . 

The discovery of the non-finiteness of the isovector spin-spin and tensor potentials due to 
class VIII in the static (M — > oo) limit points towards a further limitation of the heavy 
baryon approach when applied to the two-nucleon system. It seems that once the number of 
heavy nucleon propagators exceeds a critical value the corresponding two-loop integrals do 
not exist anymore (for M — > oo). Of course a deeper understanding of this feature would be 
very desirable. 

Finally, we consider the diagrams of class IX shown in Fig. 5. The isospin factor of the first 
and second diagram is —6 — T\-t-i while that of the third and fourth diagram is 6 — T\ ■ r 2 . The 
iterative parts of the last two diagrams are separated off in the standard way by switching 
one sign of an e-term in a single nucleon propagator. After a lengthy calculation we find for 
the imaginary part of the isoscalar central NN-amplitude Vc from class IX, 

Im ^ x, ^ = ^//^{ 1 + ^^7[|- 2arccos( - 2) ]} ■ < 21 > 

Z 2 <1 
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Fig. 5: 3ir- exchange diagrams of class IX proportional to g\. 



In order to arrive at the very short expression given in eq.(21) we have exploited the 
inherent permutational symmetry of the 37r-phase space integration. In the chiral limit m n = 
one obtains from class IX an attractive isoscalar central potential of the form 



9 A r ~ 



(4vr 2 - 60) , (22) 



which just cancels the one coming from class VIII given in eq.(18). Later when presenting 
numerical results in Table 1 we will see that this cancellation holds also for any finite pion 
mass m n . Furthermore, we obtain from class IX a non- vanishing contribution to the isovector 
central NN-amplitude Wc- For the corresponding imaginary part all appearing integrals can 
be solved analytically and one finds the following closed form expression, 

ImW#*>(.» = ^~ (4 ^ 6 )2 (3m» + 2»ml - 9^ - 4/x 3 ) . (23) 

Inserting this mass-spectrum into eq.(2) one obtains finally a repulsive isovector central po- 
tential which can be expressed in terms of a simple exponential-function multiplied by a 
polynomial, 

Wc X \r) = —^- =- (l+m T r)(4 + 5m f r + 3m 2 / 2 ) . (24) 

We have also evaluated the imaginary parts of the isoscalar and isovector spin-spin and tensor 
NN-amplitudes arising from the four diagrams shown in Fig 5. In these cases there appear 
again singular terms proportional to (1— z 2 )~ 3 ^ 2 in the integrands such that the double-integral 
Jj z 2 <i dui\dui2 diverges. Consequently, we may conclude that only the central potentials gen- 
erated by the 37r-exchange diagrams proportional to g\ (shown in Figs. 4,5) exist in the static 
limit (M — > oo) while all spin-spin and tensor potentials diverge. This is a somewhat unex- 
pected feature of the 37r-exchange NN-potential. 

In Table 1, we present numerical results for the coordinate space NN-potentials generated 
by the 37r-exchange graphs of class V, VI, VII, VIII and IX for internucleon distances 0.6 fm < 
r < 1.4 fm. We use the parameters f n = 92.4 MeV, m n = 138 MeV (average pion mass) and 
9a — g^Nf-ir/M = 1.32 employing the Goldberger-Treiman relation together with g nN = 13.4. 
The choice Qa = 1.32 is most natural in the present context since the pion-nucleon coupling is 
the relevant here and not the axial-vector coupling. One observes that the isovector spin-spin 
and tensor potentials from classes V and VI (scaling with g\) are a factor 3 to 10 larger than 
the (largest) potential W^t ( r ) found in our previous work [1] on the chiral 37r-exchange 
diagrams scaling with g\. The most striking result is that the isoscalar central potentials 
from class VIII and IX cancel each other, Vq 1 (r) + V^ IX \r) = 0, and this cancellation is 
found to happen with high numerical precision. From the double-integral representation in 
eqs. (17,21) it is not at all obvious that the relation lmV^ VIII \ijj) + lmV^ IX \ijj,) = holds, at 
least we have not yet found a simple analytical proof. Consequently, we can conclude that the 



total isoscalar central potential generated by chiral 37r-exchange vanishes identically. Only 
the exchange of an explicit cj(782)-vector meson (a strongly resonant 37r-state) can contribute 
to the (phenomenologically needed) isoscalar central repulsion. 

By inspection of Table 1 we furthermore observe that the largest chiral 37r-exchange po- 
tentials are the repulsive isoscalar spin-spin and tensor potentials V^p H \r) due to class VII 

— — ( J V\ 

and the repulsive isovector central potential Wq (r) due to class IX. These are just the cases 
in which simple analytical expressions could be given in eqs.(15,16,24). Actually, we are dis- 
proving here the claim of ref . H] that a particular diagram of class I (namely the first graph in 
Fig. 1 of ref. [jl|) would be dominant. Presumably, the strength of the potentials V^. (r) and 

(r) is still too weak to be of practical relevance. Of course, in order to confirm this an 
empirical analysis of the low-energy NN-data base as done in ref. including these dominant 
chiral 37r-exchange potentials would be very desirable. 

In summary, we have completed in this work the calculation of the static chiral 37r-exchange 
NN-potentials by evaluating all diagrams proportional to g\ and g\. As a surprise we find that 
the total isoscalar central potential vanishes identically. The individually largest chiral 37T- 
exchange NN-potentials are of isoscalar spin-spin, isovector central and isoscalar tensor type 
and for these potentials very simple analytical expressions could be given. We furthermore 
discovered that the spin-spin and tensor potentials due to class VIII and IX (consisting of 
diagrams proportional to g\) do not exist in the heavy nucleon mass limit M — ■> oo. 



r [fm] 


0.6 


0.7 


0.8 


0.9 


1.0 


1.1 


1.2 


1.3 


1.4 


W ( S V) [MeV] 


51.7 


16.5 


6.02 


2.43 


1.06 


0.503 


0.249 


0.129 


0.069 


[MeV] 


16.0 


4.99 


1.79 


0.704 


0.302 


0.141 


0.068 


0.035 


0.018 


W [ s 1] [MeV] 


19.9 


5.99 


2.06 


0.783 


0.322 


0.141 


0.065 


0.031 


0.015 


W? I] [MeV] 


-30.2 


-9.78 


-3.63 


-1.52 


-0.689 


-0.326 


-0.165 


-0.087 


-0.048 


V^ vn) [MeV] 


214.23 


68.81 


25.38 


10.40 


4.62 


2.19 


1.10 


0.575 


0.313 


% VII) [MeV] 


58.39 


19.20 


7.25 


3.05 


1.39 


0.675 


0.347 


0.186 


0.104 


V$ lll) [MeV] 


177.6 


56.2 


20.4 


8.21 


3.58 


1.68 


0.830 


0.429 


0.229 


V^ X) [MeV] 


-177.6 


-56.2 


-20.4 


-8.21 


-3.58 


-1.68 


-0.830 


-0.429 


-0.229 


W { c X) [MeV] 


148.80 


46.54 


16.75 


6.70 


2.91 


1.36 


0.665 


0.342 


0.183 



Tab.l: Numerical values of the local NN-potentials generated by the chiral 3tt -exchange graphs 
of classes V, VI, VII, VIII, IX (shown in Figs. 1, 2, 3, 4, 5) versus the nucleon distance r. The 
units of these potentials are MeV. 
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